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We present calculations of the decay constants and kaon B-parameter Bk as the first stage of RBC Collabo- 
ration's quenched numerical simulations using DBW2 gauge action and domain-wall fermions. Some of potential 
systematic errors and consistency to previous works are discussed. 



1. Introduction 

In the quantities related to kaon physics such as 
/,,-, Bk, ReA /KeA 2 and e'/e, the first three 
are constructed straightforwardly on the lattice. 
However, in actual numerical simulations, there 
are several potential sources of systematic error, 
namely explicit breaking of chiral symmetry, scal- 
ing violation, finite volume effect and quenching 
effect. Although many efforts have been made to 
decrease these errors, yet more numerical simula- 
tions should be done to attain a conclusion. 

We report our results of f^, Jk and Bk as a 
conclusion of the first stage of our longstanding 
quenched numerical simulations using domain- 
wall fermion and DBW2 gauge action. By our 
choice of the lattice action, it is expected that the 
contamination from the explicit chiral symme- 
try breaking is vanishingly small. In particular, 
we estimate the effect of operator mixing in the 
calculation of Bk through the non-perturbative 
calculation of the renormalization factors which 
weight contributions of the mixing operators. We 
also examine the scale dependence of our results 
by carrying out two kinds of numerical simula- 
tions with a similar physical lattice volume « 1.6 
fm but with different scales a" 1 w 2 GeV and 3 
GeV. After comparing our results with those from 
previous works, we present our continuum results 
of B K . 

For other kind of RBC's numerical simulations 
with the dynamical quark, see ref. £Q. 



*We thank RIKEN, BNL and the U.S. DOE for providing 
the facilities essential for the completion of this work. 



Table 1 

Simulation parameters, statistics and results of 



a 1 , m rcs a, 


Za and m s a/2. 







1.22 


1.04 


size 


24 3 x 48 


16 3 x 32 


M 5 a 


1.65 


1.8 


L s 


10 


16 


rrifa 


0.008 - 0.040 


0.01 - 0.05 




in step of 0.008 


in steps of 0.01 


^configs. 






(Zb k ) 


53 


50 


(others) 


106 


202 


basic results 




a- 1 


2.914(54) GeV 


1.982(21) GeV 


m TCS a 


9.73(4) • 10~ 5 


1.85(12) • 10~ 5 


Z A 


0.88813(19) 


0.84019(17) 



2. Numerical Simulations 

In Tabled we enumerate two sets of our sim- 
ulation parameters, statistics and results of the 
basic quantities such as lattice scale from the in- 
put m p = 770 MeV, the residual quark mass m res 
and the renormalization factor of axial vector Za- 
Results of the basic quantities for j3 ~ 1.04 are 
quoted from ref. [5], in particular. For the finer 
lattice with (3 = 1.22, it is known that the topo- 
logical charge Qtop changes very slowly in the or- 
dinary Markov chain [2]. To avoid incorrect dis- 
tribution of Qtop; we generated configurations as 
described in ref. 0. As a result, we obtained a 
reasonable distribution of the topological charge: 
(<3to P ) = 0.38(29). We employ the averaged 
quark propagator over those with periodic and 
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Figure 1. Results of / w (filled symbols) and 
[GeV] (open symbols) versus a 2 [GeV -2 ]. 
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Figure 2. Zb k and RG1 value Zb k versus p 2 att 
for DBW2 P = 1.22 (left) and = 1.04 (right). 
In each panel, Zb k (circles), RGI value (squares) 
and its linear extrapolation are shown. 



anti-periodic boundary conditions in the tempo- 
ral direction for the calculation of / W! fx and Br. 
Therefore, temporal lattice size may be treated as 
96 and 48 for /3 = 1.22 and 1.04. 

3. Decay Constants 

As a combination of fit parameters of two point 
correlation functions of mesons, we calculate de- 
cay constants of the pseudo-scalar from the com- 
bination of the fit parameters: 



/raps."!/ APP 

y 2 -^ww 

where Ap* p and are the amplitudes of 

the correlation functions (Al° int {t)P wM (0)^ and 

(P wall (t)F wall (0)). mps is determined from the 
simultaneous fit of these two correlation functions 
thus common to both. Using the linear fit with 
fps — co + cimpg, we extract bare values of f n at 
mps — rn^ and Jk at mps = ttik- Multiplying 
these values by Za in Tabled we show our phys- 
ical results in Figure ^ As well as for j3 = 1.22 
(circles) and 1.04 (squares), results for (3 = 0.87 
(diamonds) from ref. [2] are plotted in this figure. 
By the linear fit using these three data, our results 
are U = 140.1(3.3) MeV and f K =153.9(2.6) 
MeV in the continuum limit. They are both in- 
consistent with the experimental values. On the 
other hand, our result Jk/ f-K = 1.098(13) in the 
continuum limit is roughly consistent with the es- 



timation from the quenched chiral perturbation 
theory 0|. 

4. Kaon B-parameter Bx 

We compute Bps i.e. Bk on the lattice, by 
extracting a plateau of the ratio of the correlation 
functions: 

(0|P wall (^)Q A g = 2(t)P wall (t O )|0) 
| (0 IPwaii^jA^)! 0) (0 \A 4 (t)P^(t )\ 0) ' U 

where the locations of two wall sources P wa11 are 
(t ,t' ) = (7, 41) and (5, 27) for f3 = 1.22 and 1.04. 
We chose 19 < t < 29 and 14 < t < 17 for the fit 
ranges in each case. 

To determine the renromalization factor Zb k , 
we employ the non-perturbative calculation 
which was pioneered by ref. Including the 
possibility of the contribution from the mixing 
operators, the renormalization condition is writ- 
ten by 

z- 2 ^rg latt (p latt ) = rg tree , 0) 

where p^ Iatt a nd r^'' tree are the amputated four 
point vertices on the lattice and in the tree level 
for the four-quark operator Oj with certain chi- 
ralities i = VV ± AA, SS =F PP and TT. By 
solving (0, we observe that the renormalization 
factors Zij for the mixing operators are at most 
- 0.1% of the factor for VV + AA, which cor- 
responds to Qas=2- These small Zy's suppress 
the effect of the operator mixing to be less than 
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the statistical error of Bps in spite of the fact the 
B-parameters of the mixing operators are a few 
dozens times larger than Bps 

In Figure [3 we show the results Zg I ^ MOM = 
Zvv+aa,vv+aaI Z\ hi RI/MOM scheme and the 
renormalization group independent (RGI) value 
(piatt)Z RI/MOM b la tt), where 



w 



RI/MOM 



the factor w^a/MOM was calculated in ref . [7| . Ex- 
trapolating the data for pi a tt > 1 linearly, we 
quote Zb k at piatt = 0. Renormalization fac- 
tor in the MS, NDR scheme are calculated as 
MS ( M = 2GeV)i Bjf . 
We plot renormalized value B ps 



Bk 



2 GeV)Bps as a function of mp s [GeV ] in Fig- 



ure|3 where circles are from (3 — 1.22 and squares 
from 1.04. Using the chiral expansion [5| 
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and our results of fp s(mps — 0) for the co- 
efficient of the chiral log, we obtain the solid 
and dashed curves for (3 = 1.22 and 1.04 with 
^ 2 /dof = 0.80 and 0.08, respectively. Result of 
B^ s (p = 2 GeV) for each (3 from the interpo- 
lation to mps — rriK are shown by the filled 
symbols in Figure 0| In the same figure, we plot 
results of CP-PACS 0| using Iwasaki gauge ac- 
tion with two set of similar parameters as ours, 
i.e. {a- 1 = 2.87 GeV, 24 3 x 60, L s = 16) 
and {a- 1 = 1.88 GeV, 16 3 x 40, L s = 16) 
with the open squares. On the other hand, the 
open circle corresponds to the previous result of 
RBC ^0] using Wilson gauge action with (a" 1 = 
1.922(40) GeV, 16 3 x 32, L s = 16). For a" 1 « 
2 GeV, results from these three kinds of calcula- 
tions distribute over the width of more than 10%. 
On the other hand, for a -1 3 GeV, results from 
both Collaboration agree with each other. In the 
continuum limit taken by the naive extrapolation 
of our two results, we extract B^ s (^, — 2 GeV) = 
0.569(21) and B K = 0.762(28). 
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Figure 3. Bpg vs. m 2 [GeV^] with fit curves. 
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